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Abstract 

We show that complete cotangent lifts of vector fields, their decomposi- 
tion into vertical representative and holonomic part provide a geometrical 
framework underlying Eulerian equations of continuum mechanics. We 
discuss Euler equations for ideal incompressible fluid and Vlasov equa- 
tions of plasma dynamics in connection with the lifts of divergence-free 
and Hamiltonian vector fields, respectively. As a further application, we 
obtain kinetic equations of particles moving with the flow of contact vec- 
tor fields both from Lie-Poisson reductions and with the techniques of 
present framework. 
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1 Jets 

Let (£, 7T, Ad) be a smooth bundle with coordinates (x a ; 1 < a < dim (Ad) = to) 
on the base manifold Ad and (x a , u x \ 1 < A < rank (tt) = k) on the total mani- 
fold £ . The vertical bundle associated with tt is 

Vtt = kerTir = {£ G T£ : Ttt (£) = 0} (1) 

and this is a vector subbundle of the tangent bundle T£. Here Ttt denotes the 
tangent mapping of the projection tt. Two sections cf>,ip G & (tt) of the bundle 
tt at a point x G Ad are called equivalent if their tangent mappings are equal at 
that point, that is, T x <f> = T x tp. Given a point x, an equivalence class containing 
a section <j) is denoted by j x 4>- The first order jet manifold 

J 1 tt = {jl(j> : x G M and <j> G © (tt)} (2) 

associated with (£, tt, Ad) is the set of equivalence classes at every point x G Ad 
with induced coordinates 



(x a ,u\u x a ) : J 1 TT^R m + k+mk -.j 1 ^^ x a , U A (0(x)),^ 



(3) 



We have fibrations ttq : J 1 tt £ : j^cj) — > <f> (x) and tt\ : J 1 tt Ad : j x (j> —> x of 
J x tt on £ and Ad, respectively [5], [15]. 
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Given a differentiable map p : JV — > A4 from a manifold J\f to the base 
manifold Ai, the pull-back bundle of it by p is the triple (p* £ , p* n , Af) where 

p*£ = N x M £ = {(n, e) e N x £ : tt (e) = p (n)} (4) 

is the Whitney product and, p*7r = pr\ is the projection to the first factor [SJ. 
Consider the pull back bundle 

(ttS (T£) = Ax £ T£, 7t* t £ = pn, J 1 *) 

of (T£,T£,£) by the projection ttq : J 1- 7r — > £, where Tg is the tangent bundle 
projection. A section of itqT£ is called a generalized vector field of order one 
[IS], [IB]. One may regard a section of ttqT 5 as a map from J 1 ^ to T£. We 
require that generalized vector fields are projectable [BJ. 
In coordinates, a generalized vector field is 



and its first order prolongation pr 1 ^ is 



(5) 



where D x a is the total derivative operator with respect to x a and, u x a (j x <t>) = 
is an element of the second order jet bundle. Lie bracket of two 
first order generalized vector fields £ and r\ is the unique first order generalized 
vector field 

% V] pr0 = (Pr 1 ^ Vr\ (O) ^ + {pr'Z tf) W \ {t)) ^ (7) 

If £ and 11 are two vector fields on £, then [ , 1_„ reduces to the Jacobi-Lic 

a I I L ' 1 JJTO 

bracket of vector fields [13] ■ 



2 Lifts 

Consider a vector field X e X(A4) on A4, and let be a section of 7r. The 
holonomic lift of X (x) e X^A^ by is 

(X (x))) e ttJ TO - x £ T£. (8) 

In coordinates, if A = X a (x) d/dx a , then 

* M - + - + **■ (9) 

Define the holonomic part of a projectable vector field £ G X (£) as the holonomic 
lift of its push forward by 7T, that is 

Hi = (n^) h ° l ■ (10) 
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H£ is a generalized vector field of order one. Define a connection (1; 1) tensor 



satisfying H£ — Lj£. Then, the vertical (or evolutionary) representative 

vz = a - rj (o = (r - r^) ^ (12) 

of £ is vertical valued generalized vector field of order one [13] , [H] , [IS] ■ 

Proposition 1 Holonomic lift is a Lie algebra isomorphism from the space of 
projectable vector fields in X(£) into J l it x$ T£. 

Proof. We consider two projectable vector fields £ and r\ on £ . A straight 
forward calculation gives 

[rj (0 , Tj m pro = k hol ,V h ° l ] = K, rj\ ho1 = T 3 K, rj] (1 3) 

L J pro 

where [ , ] pro is i/ie Lie bracket for generalized vector fields in Eq.^Tty. ■ 

On the other hand, the generalized bracket of vertical representatives satisfies 

[V^Vv} pro = V[^v] pro + < B^,v), (14) 
where 25 is a vertical- vector valued two-form 



[v hol ,V(\ - e o \Vr, . (15) 

^ . . WTO 



There is, however, a class of vector fields, defined again by lifts, for which the 
vertical representative becomes a Lie algebra isomorphism. Let (p t : M. — > M. 
be the flow of X on M.. Cotangent lift of <p t is a one-parameter group of 
diffeomorphism tp" on T*Ai satisfying 

TT M O Iff* =(fi t OTT M (16) 

where ttm is the natural projection of T*M to A4. The cotangent lift of 
the inverse flow T*ip_ t satisfies the argument in Eq. (|16[) . Infinitesimal gen- 
erator X c * : T*M -> TT*M of the flow Lpf* is called complete cotangent lift 
of X. X c * is a Hamiltonian vector field on the canonical symplectic manifold 
(T*M, Qt*m = —dOT'M) for the Hamiltonian function P (X) = i X "9r*M M- 
The infinitesimal version 

Tttm ° X c * = X o it M ■ 

of Eq. (fT6"|) gives the relation between X and X c * with T7r» being the tangent 
mapping of ttm- The complete cotangent lift mapping c * : X(M) — > X(T*A4) 
taking X to X c * is a Lie algebra isomorphism into [8].|17j 

[X C *,Y C *] = [X,Y] C * , VX,Y eX{M). (17) 
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In Darboux's coordinates (x a ,y,b) on T*A4, the complete cotangent lift of 
X = X a (x) d/dx a on M is 

d dX b d 

X " = XViX) = ~ Vb ^Wa (18) 

with the Hamiltonian function being V (X) (x, y) = ybX b (x). We decompose 
the complete cotangent lifts into vertical representative and holonomic part 

where the connection in Eq. (|llj) has the particular form 

r -*-«G* + £<s)- (20) 

Proposition 2 TTie mapping V c * : X(M) -> X(T*X) : X -> ^X c * is a Lie 
algebra isomorphism into. 

Proof. The vector valued two form *B in Eq. H15\) vanishes for the complete 
cotangent lifts, that is, 03 (X C *,Y C *) = for all X, Y e X(M) , therefore one 
has V [X C *,Y C *] = [VX C *, VY c *} pro and the result 

V[X,Y] C * = [VX",VY c *] pro . (21) 

follows from Eq.{l r fy. ■ 

The last object we consider in this section is the vertical lift of one forms. 
Take the cotangent lift T*n M : T*M -> T*T*M of the projection n M : 
T*M -> M and recall the isomorphism ^ T . M : T*T*M -)• TT*M associ- 
ated with the symplectic two- form £It*m on T*M. Define the Euler vector 
field 

X E : T*M -> TT*M : z -> n*,.^ o T*tt m (z) (22) 
which is vertical, that is, image (X E ) C ker (Tttm)- Indeed, 



{z,Tn M o X E (z)} = (T*n M (z) ,n\,, M oT*n M (z)^ 

= !J T . M (n M (z),ri M (z))=o, (23) 

Vz G T*A4, where we used the skew-symmetry of VLt'm- X e is the unique 
vector field satisfying the following equalities 

ix E ^T*M — Ot'Mi £x e ^T'M — —Qt»m, £-x e 6t'M — ~0t*m, (24) 

where ix E and Cx E are the interior product and the Lie derivative operators 
[7]- Let aeA 1 (M) be a one-form on M. The vertical lift 

a v = X E o a o tt m : T*.M -> TT*X (25) 
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of the one- from a is a vertical vector field on T*A4. The Jacobi-Lie bracket of 
a complete cotangent lift and a vertical lift is a vertical lift 

lX c *,a v ] = (C x a) v (26) 

for X G X(M) and a G A 1 (A4) [17]. In coordinates (z a ,y 6 ) of T*M, the 
Euler vector field is Xe — —y a d/dy a and the vertical lift of the one-form a = 
a a (x) dx a becomes a v = —a a (x) d/dy a - 

3 Dynamics 

Assume that a continuum initially rests in M, and the group Diff(M) of 
diffeomorphisms acts on left by evaluation on M. 

Diff(M)xM^M:(<p,x.)^ip(x.) (27) 

to produce the motion of particles. The right action of Diff(Ai) commutes 
with the particle motion and constitutes an infinite dimensional symmetry group 
of the kinematical description. This is the particle relabelling symmetry [2]. 
An element of the tangent space T v Diff(M) at ip G Diff(M) is a map 
V v : M. — > TA4 called the material velocity field and satisfies r^vi °V v =(p. In 
particular, the tangent space T it i M Dif f (A4) at the identity i<1m & Diff(A4) 
is the space X (A4) of smooth vector fields on M. The Lie algebra of Diff(M) 
is X(Ai) with minus the Jacobi-Lie bracket of vector fields [8]. 

The dual space X* (M) ~ A 1 (M) <8> Den(M) of the Lie algebra is the 
space of one- form densities on M. The pairing between a <g> dfi G X* (A4) and 
A G X (A4) is given by 

{a®d[x,X)= [ (a(x),A(x))d^(x). (28) 

J M 

The pairing inside the integral is the natural pairing of finite dimensional spaces 
T X M and T*M. The coadjoint action is 

ad* x : X* (M) -> X* (A4) 

: a ® — )■ £x (a ® d/i) = (£xo + (diVd^X) a) <8> dyu. (29) 

VA G X (A^) and hence the Lie-Poisson equations on X* (M.) are 

& = -C x a - (div d ^X) a, (30) 

where divd^X denotes the divergence of the vector field X with respect to the 
volume form dfi. 

In terms of vertical lifts, the dynamics in Eg. (1301) is generated by the vector 
field {Cxot + (divdfiX) af . For the divergence free vector fields, if a = y a dx a , 
then the Lie-Poisson equations are generated by 



{C x (y a dx a )) v = VX°* (x a ,y a ) 



(31) 
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3.1 Ideal incompressible fluid 

For an ideal incompressible fluid in a bounded compact region QcK 3 the con- 
figuration space is the group Dif f vo i (<2) of volume preserving diffeomorphisms 
on Q. The Lie algebra %div (Q) of Diff vol (Q) is the algebra of divergence free 
vector fields parallel to the boundary of Q and, the dual space X^ iv (Q) is the 
space 

(S) = {[T] ® d 3 q G (A 1 (Q)/dJ-(Q)) ® £ e n(Q)}, (32) 

of one-form modulo exact one-form densities on Q. Here, [T] = {T + dp : p G J-(Q)} 
denotes the equivalence class containing T and the volume three form d 3 q is 
the Euclidean volume on R 3 [2] , [IT] . 

Let (x a , T;,) be induced coordinates and X = X a d/dx a be a divergence free 
vector field. The complete cotangent lift of X is 

X c * = X a J^ - T b (dX h /da" s 



and its vertical representative becomes 

Equations of motion for the dynamics generated by VX C * are 

3[T1 



dt ="£xPl. (34) 

For a generic element T + dp G [T] , Eq. (|34f becomes Euler's equations for ideal 
fluid, that is dT /dt + Cx^ = dp. If the dual space X^L (Q) is identified with 
exact two forms by [T] — > c?T = weA 2 (Q), then Eg. becomes the Euler's 
equation in vorticity form dw/dt + Lx^> = 0. 



3.2 Collisionless plasma 

We take Ai to be cotangent bundle T* Q of Q C K 3 in which the plasma particles 
move. The configuration space of collisionless nonrelativistic plasma is the group 

Diff can (T*Q) = { V G T*Q : ^ T . s = n T « Q } (35) 

of all canonical diffeomorphisms where Qt*q is the canonical symplectic two 
form on T*Q We assume that, the Lie algebra of Diff can (T*Q) 

is the space of globally Hamiltonian vector fields X^ am (T* Q) with minus the 
Jacobi-Lie bracket so that the equations 

[Xh,X f ] JL = -X{h,f}n T . a ( 36 ) 

describe a Lie algebra isomorphism 

h^X h :(T(T*Q),{, }n T . a ) ^(Xham(T*Q),-[, } JL ) , (37) 

between %ham (T* Q) and the space of smooth functions J- (T* Q) modulo con- 
stants endowed with the (nondegenerate) canonical Poisson bracket { , }n T , Q - 
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Proposition 3 The dual space of the Lie algebra Xham {T*Q) of Hamiltonian 
vector fields is 

*Lm - Wrf ® d V 6 A X (T*Q) ® Den(T*Q) : dwn r , e n? d ^ 0}. (38) 

With this definition of the dual space the -^-pairing of the Lie algebra and 
its dual becomes nondegenerate provided we take the volume form to be the 
symplectic one d/i = Oy» g in 

/ (X h (z) , n td (a)) d M (z) = - / (dh,U i id )dn=- [ i ni {dh)dfi 
Jt"Q Jt*q * ' jt*q * d 



— / dhAi n } d/i = hdintdfi 

Jt*Q id JT'Q 

I hdivn T * Q Ill d dti, (39) 



where we use the musical isomorphism q : Ily —> H id induced from the sym- 
plectic two- form Qt* q and apply integration by parts [5J internet supplement] . 
The dual of the Lie algebra isomorphism in Eq. (l3"Tl) is 

U id (z)^divn TrQ Ul( Z ) (40) 

and it is a momentum map. In Darboux's coordinates z = (q l ,Pi) on T*Q, 
we have Ot*q = dq l A dpi and we take = IL (z) dq l + IF (z) dpi. Then, the 
momentum map 



tt\ a- tt(j / \ dw (z) dn. t (z) 

/ (z) - dn^n?, (z) = ^— (41) 



defines the plasma density function. 

In the induced coordinates (q l ,pj; Hi, IP) on T*T*Q, consider the Hamil- 
tonian function h = (l/2m) 5 lJ piPj + e</> (q) which is the energy of a charged 
particle on Q |]. The corresponding Hamiltonian vector field is 

v i \ 1 Hi ^ d(j> d 

X h (z = — ^Pi^-r " e JTi~^ — • 42 
m oq 3 oq % dpi 

The complete cotangent lift of Xh and its decomposition into vertical represen- 
tative and holonomic part are 

VK = (^^-Xm))~^ + X h ^))^, (43) 
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where Xh (IT) denotes the action of Xh on LL. Since, Hamiltonian vector fields 
are divergence free, the Lie-Poisson equations 

if = -x h (it) - — (f J n 7 - (44) 

are generated solely by VX£*. These are Vlasov equations in the momentum 
variables [I]. For the density formulation, we make back-substitution of the 
plasma density function / (z) = diva T , Q Tl\ d into Eqs. (|44|) and obtain the 
Vlasov equation 

df | Ppj df Q d<pdf =Q 
dt m dqi dq l dpi 

3.3 Contact flows in 3D 

Let M. be a three dimensional manifold with a contact one form a £ A 1 (Ai) 
satisfying da A a ^ 0. A contact form determines a contact structure which, 
locally is the kernel of the contact form a. A diffeomorphism on A4 is called 
a contact diffeomorphism if it preserves the contact structure. We denote the 
group of contact diffeomorphisms by Diff con (M). A vector field on a contact 
manifold (M,a) is called a contact vector field if it generates a one-parameter 
group of contact diffeomorphisms [Tj.|12j. 

In Darboux's coordinates (x,y,z) on A4, we take the contact form to be 
a = xdy+dz. For a real valued function K — K (x, y, z)onM, there corresponds 
a contact vector field 

X -(—- x —\—-——+(-K + x—\— (46) 
\ dy dz J dx dx dy \ dx ) dz 1 

on A4 satisfying the identities 

ix K & = —K and ix K dcr = dK — [i^dK) a, (47) 

where R a = d/dz is the Reeb vector field of a. R a is the unique vector field 
satisfying iu a a = 1 and i^da — 0. The divergence div^^X^ of X^ with respect 
to the volume form dfi = da A a can be computed to be div^XK = —2R a K. 

Contact Poisson (or Lagrange) bracket of two smooth functions on A4 is 
defined by 

_ dL dK dL 3K | dK / . x 8L\ _dLf R _ x dK\ 
' c dx dy dy dx dz \ dx J dz \ dx J 1 

\/L,K £ T (M) . The identity [Xk,Xl] jl = —X{ KL } establishes an isomor- 
phism between Lie algebras (X con {M) , — [ , ] JL ) and (J 7 (M.) , { , } c ). Follow- 
ing result gives a precise definition of the linear algebraic dual of X con {M)- 



9 



Proposition 4 The dual space of the algebra X con (M) of contact vector fields 



Kon {M) = {ac^dfieA 1 (M) (g> Den (M) : da A a - 2a A da ^ 0} (49) 
where a is the contact form on M. and dfi — da A a . 

Proof. Proof. This follows from the requirement that the pairing between 
X con and X* on (-A4) be nondegenerate. We compute 

/ (a,Xi()dn = / a/\ix K da/\a+ I (ix K a)aAda 
JM JM JM 



a A (dK - (i Rtr dK) a) A a 



M 



Ka A da 



M 



I K (da A a -2a A da) 
Jm 



(50) 



where we use the identities in Eq. J^7|) at the second step 



A geometric definition of density of contact particles can be achieved by 
considering the Lie algebra isomorphism T (M) — ► X CO n (M) ■ K — > the 
dual of which is a momentum map 



Kon (M) 0^) ■ a ^ da Aa - 2a Ada. 

and defines a real valued function L on 

Lder A c = da A cr — 2a A dcr. 



(51) 



(52) 



In coordinates, let a = a x dx + a y dy + a z dz G 5E* „ (-A4) and recall cr = xdy + dz. 
Then, 



L(x,y,z) 



da x da y 
dy dx 



da z da x 
x- 



dx dz 

The dual space X* on (TV) admits the Lie-Poisson bracket 

~5F) 58. 

5a ' 5a JL 



2a z . 



(53) 



{f),&}(a) = - f 

JM 



a, 



dp = - (a,[X H ,X K ] JL )dn, (54) 

JM 



where ft, Si <E J 7 (X* con (M)) and ^/(5a = X H , 5Si/5a = X K e X con (A4). The 
Hamiltonian operator J^p (a) associated to the Lie-Poisson bracket in Eg. (jM]) 
is defined by 

{Sj,Si}(a) = - [ (X H ,J LP (a)X K )dfx (55) 

JM 

and a direct computation gives 

Proposition 5 The Hamiltonian differential operator associated to the Lie- 
Poisson bracket in Eq. \5J$ is 



Jlp (a) 



(3d d d d d \ 

1 a. x — + — ■ a x Oi y — + — -a x a z — + — -a x \ 

ax ox ox ay ox Oz 

d d do d d 

a xs~ + ■ ay ot y — + — • a y a z — + — ■ a y 

Oy Ox oy oy Oy Oz 

O d O d 

\ct x — + — -a z a y — + — -a z a z — + — ■ a z i 

\ oz Ox Oz oy Oz Oz / 



(56) 
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where d/dx ■ a y — a y d/dx + da y /dx. Assuming 5M./5a = Xk, the Lie-Poisson 
equations on 3£* on (Ad) are 

a = Jlp (a) X K = -ad* XK a = -Cx K a - (div^X K ) a. (57) 

The Lie-Poisson bracket on the dual space Den (Ad) of J 7 (Ad), as defined 
by Eq.flSTJl, is 

{M( L )=J L {^^} d»= J M L{H,K} c d^, (58) 

where Sfi/SL = H, 5&/5L = K € T (Ad) and d[i = da A a. 

Proposition 6 The Hamiltonian operator Jlp (L) for the Lie Poisson bracket 
in Eg. J5^|) is 

dL\ d 



J LP (L) = X L+ [AL + ^) Tz , (59) 

and the Lie-Poisson equation on Den (Ad) becomes 

L = -{L, K} c - 2div^ (X K ) L. (60) 

Proof. The verification of the Hamiltonian operator in Eg. i59\) is a straight- 
forward calculation which follows directly from the definition of the Lie-Poisson 
bracket in Eg. J5^| ). To obtain the Lie-Poisson equation we compute the coadjoint 
action negative of which is the reguired eguation. By definition 

(ad* K L, H) = (L,ad K H) = (L,{K,H} c ) 

= - J^L{H,K} c dfi=- J^L^X K (H) + ^H^d^ 

= J (x K (L) + div dt , (X K ) L - Hdfi 

= J^^{L,K} c -^L + div dfl (X K )L-^pjHd f x 

= f ({L,K} c + 2div dfi (X K )L)Hdfx, (61) 

where we use integration by parts at the third step and the identities 

{H, K} c = X K (H) + = -X H (K) - ^-K (62) 

oz oz 

at the second and fourth steps. ■ 

The equation of motion L = —ad* K L is the kinetic equation of contact par- 
ticles in density formulation. 

Proposition 7 The Hamiltonian differential operators Jlp (a) in Eg. A56\) and 
Jlp (L) in Eg. H59}) are related by 

HJ LP (L)K = -X H J LP (a)X K (moddiv). (63) 
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We now obtain dynamics of contact particles by the methods of previous 
sections. Let (A4, a) be a contact manifold and consider the contact vector field 
Xk in Eq. pij)) . Its complete cotangent lift is 

v „ v f^dK dK\ d T J dK 9 dK 8 \ raA . 

x -= XK+ { r ^ + ^)^ + ^ + ^{-^^- y + ^o^) (64) 

where we use the following abbreviations 



9 \ d d d d 

;— , V = a y - a x - xa z — , $ = xa x — 

ox ) ox dy ox dz 



and the induced coordinates (x,y,z,a x ,a v ,a z ) on T*J\f. X^ is a canonically 
Hamiltonian vector field. The vertical representative VX^ of X^ is 



oz ox J da x \ oy J oa y 

dK \ d 

(H$)^-I K (a z ) (66) 
oz J oa z 

with Xk {oi x ) denoting the action of Xk on a x . To obtain the equations of 
motion for the momentum variables, one needs to add the divergence term, 
that is, 

a = VXj? (a) - (div dtl X K ) a. (67) 

It can be checked that Eq. tjSTf and Eq. (j57l) are equal. In coordinates, the system 
of equations in Eq. (|57| takes the form 

dK dK . . dK 

a x = T— + ^—-X K {a x ) + 2—a x 
oz ox oz 



dK dK 
= (t> + 1 *)^- X K(a v )+2—a y 

dK , , dK 



($ + - X K (a z ) + 2— a z . (68) 

dz dz 



Substituting L in Eq.(|53f to the system of Eqs. ([68|) we obtain the evolution of 
the density of contact particles as given by Eq. ([60| . 
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